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SUMMARY 

A method for predicting Isothermal plastic fatigue life of a composite 
lamina Is presented In which both fibers and matrix are Isotropic materials. 

In general, the fatigue resistances of the matrix, fibers, and Interfacial 
material must be known In order to predict composite fatigue life. Composite 
fatigue life In this paper Is predicted using only the matrix fatigue resist- 
ance due to Inelasticity micromechanisms. The effect of the fiber orientation 
on loading direction Is accounted for while predicting composite life. The 
application is currently limited to Isothermal cases where the Internal ther- 
mal stresses that might arise from thermal strain mismatch between fibers and 
matrix are negligible. The theory Is formulated to predict the fatigue life 
of a composite lamina under either load or strain control. It Is applied cur- 
rently to predict the life of tungsten-copper composite lamina at 260 °C under 
tension-tension load control. The calculated life of the lamina Is In good 
agreement with available composite low cycle fatigue data. 


INTRODUCTION 

Fatigue failure of a metal matrix composite (MMC) Is a complex process. 
Failure modes can depend on the applied load, the properties of the matrix, 
fibers and Interface, fiber volume fraction and orientation, and the service 
temperature and environment. Depending on these factors, the active failure 
mode can be matrix dominated, fiber dominated, fiber/matrix Interfacial fail- 
ure, or self-similar fatigue damage (ref. 1). A fatigue life prediction method 
has to consider the most active mode of failure to obtain a good estimate for 
life. 


In this work a tungsten reinforced copper composite containing unidirec- 
tional continuous fibers was studied. For specimens whose fibers were oriented 
parallel to the load axis, the failure mode was matrix dominated. Therefore, 
for the ensuing analysis this mode of failure was assumed to dominate for all 
angles of fiber orientations. 

Further, It Is assumed, for simplicity, that the composite fatigue failure 
can be considered as a sequence of two events where complete fatigue failure of 
the matrix is followed by Immediate complete failure of the fibers. Matrix 
cracking introduces additional cyclic axial and shear stresses to the fibers. 
The localized nature of these stresses severely reduces the fiber residual 
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life; therefore, Its contribution to composite fatigue life can be neglected. 
The analysis in this paper is limited to composites with a strong fiber-to- 
matrix bond and no interfacial phase. In tungsten-copper composites the bond- 
ing is excellent, and the two constituents are mutually Insoluble. 


CONSTITUTIVE RELATIONSHIP 

The state variable constitutive relationship used in this paper is 
derived in (ref. 2) and is represented by the following two equations for 
inelastic strain rate and the rate of evolution of the state variable X, 
referred to herein as "resistance to flow," 

ejj = AO (J 2 /X 2 ) r s.. (1) 

x = F ( X) J £ r+1 = ± exp(CX+D) J 2 r+1 (2) 


where e j s is the inelastic strain tensor and where AO, C , D and r are 
material constants. Values of the constants at 260 °C are given in table I. 

J 2 is the second invariant of the deviatoric stress tensor Sjj. The sign of 
the rate of evolution is positive if the material hardens and negative if it 
softens or Is damaged. The superscript i refers to inelastic strain, and 
the dot represents the first derivative with respect to time. The function 
F(X) is a material state function that controls response and determines life. 
The state variable Increases for hardening materials and then achieves a stabi- 
lized or shake down value, X s , after which the material resistance decreases 
until catastrophic failure occurs, (fig. 1). It follows for the continuously 
applied cycles that the integration of 02 r+1 over the applied cycle increases 
with the increase of X s and decreases with its decrease. This integration 
is called the loading function L shown as a function of the number of applied 
cycles. The value of this integration is assumed constant and equal to its 
value at shake down Lc. Using this assumption and integrating equation 
(2) over the entire life, it follows that the relationship between the number 
of cycles to failure, Nf, under continuously applied cycles is related to the 
initial resistance X 0 and shake down resistances as follows: 

[2exp(-CX -D) - exp(-CX -D)-exp(-D)] ,,, 


where the resistance at failure equals zero. Equation (3) is used to predict 
fatigue life whenever L s , X s and X 0 are known. 


SHAKE DOWN RESISTANCE 

The evaluation of the shake down resistance X s of the composite consti- 
tuents requires an experimental relationship between this resistance and the 
applied stress under load control or the applied strain under strain control. 
The life prediction method used In this paper requires the evaluation of the 
shake down resistance Xc and load function L s . These values are employed 
in equation (3) to calculate the number of cycles to failure. The shake down 
resistance is calculated iteratively. A value for the shake down resistance 
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Is first assumed and used in the constitutive relationship to calculate the 
strain or stress response under load or strain control, respectively. The 
calculated response is applied to the experimental relationship between the 
resistance X s and the applied stress or strain to update X s . The updated 
resistance is then used to update the calculated response. This process con- 
tinues until the change In the calculated shake down resistance over an Itera- 
tion is negligible. The load function for a constituent is calculated from 
its constitutive relations taking the resistance constant and equal to its 
value at shake down. 


LOCAL STRESSES AND STRAINS IN COMPOSITE LAMINA 
Strain Control 

In this section the local stresses and strains In the constituents of a 
composite lamina with perfectly bonded continuous fibers are considered. An 
element of the composite lamina Is shown In figure 2(a). For convenience this 
element Is assumed to be composed of uniformly spaced square fibers and matrix 
strips and is subjected to an In-plane stress that Is produced by a biaxial 
strain control. It is required to evaluate the stress and strain at any point 
in the lamina. The evaluation begins by introducing a continuum displacement 
field, which represents the longitudinal displacements along the centerline 
of fibers/and matrix strips. The continuum longitudinal strain obtained by 
differentiating the displacement field Is assumed equal to that of fibers and 
matrix at any time. Therefore, the continuum longitudinal strain rate equals 
the fiber and matrix longitudinal strain rate. The fiber and matrix shear and 
strain In the lateral direction of the fibers are not equal to the continuum 
strain obtained by differentiating the displacement field. The equality of 
such strains contradicts the equilibrium at fiber surface due to the differ- 
ence between fiber and matrix mechanical properties. The shear and lateral 
continuum strain are smaller than those of the matrix and larger than those of 
fibers. Let us assume a linear displacement field, and consider the problem 
shown In figure 2(c) where the composite is subjected to biaxial strain rates. 
It can be shown that the continuum longitudinal, lateral, and shear strain 
rates in axes parallel and perpendicular to the fibers are given by 

• * 2 • 2 

®11c * 0 X C0S 6 + ®y s1n 6 

®22c = V in 9 + ®y C0S 0 (4) 

®12c “ ( ®x " ®y J sin 9 cos 9 

where e x and ey are the applied strain rates In the x and y directions, 
respectively; where efj c , 622c and «12c- are the continuum axial strains 
measured along the fiber direction, perpendicular to the fiber direction, and 
the shear strain In the plane of these two directions respectively; and where 
9 Is the angle between the x axis and fiber direction. The matrix lateral 
and shear strain rates are obtained from the continuum strain rates by assum- 
ing that the displacement along the centerlines of the fibers and the matrix 
strips Is equal to those obtained from the displacement field. It follows 
from figure 2(a) that the longitudinal, lateral, and shear total strain rates 
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in the matrix strip are represented in terms of the continuum strain rate and 
fiber lateral and shear strain rates, as follows: 


e 11m = «11c 


e 

22m 
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22c 


(1 
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22f 


1 2m 
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V) - 


e 12 fv 


(5) 


where v is the ratio of volume fractions, equal to Vf/v m (vf is the volume 
fraction of the fibers and v m is the volume fraction of the matrix) and 
where the matrix strains eii m , 622m. and e-j2m and the fiber strains 
e22f> and ®12f <^9- 2(b)) are measured In directions parallel to the corre- 
sponding continuum strain components. The lateral and shear fiber strains 
have to satisfy the lateral and shear equilibrium on the fiber surfaces. For 
nonlinear constitutive fiber and matrix response, these stresses are obtained 
by integrating the equations of stress rates simultaneously with the constitu- 
ents constitutive relationships. The stress tensor in both fibers and matrix 
is derived from their elastic strain tensor. Therefore, the stress rate ten- 
sor of fibers and matrix are derived from the corresponding elastic strain 
rate tensor, which is the difference between the total and Inelastic strain 
rate tensors. It can be shown that the components of the stress rate tensors 
of the matrix and the fibers are represented by the following differential 
equations : 
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whers K m * Efp/(1 - u^) and, Kf = Ef/((1— Uf ) st n d Ejjj, G^, u^, Ef, Gf and Uf 
are the modulus, shear modulus, and polsson ratio for the matrix and fibers, 
respectively. The superscripts, e and i, refer to elastic and inelastic 
strain rates, respectively. The matrix stress components Sii m , S22m and 
Si 2 m and the fibers stress components Suf, S 22 f and Si 2 f are parallel to 
the matrix and fiber strain components having the same indices. The rate of 
elastic lateral and shear strain in the fibers is obtained from the lateral and 
shear equilibrium between the matrix and fibers and can be written as follows: 

•fcf - [A2/K f ♦ v V ] - [ Uf K f /(K f , »K m >Ke 11c - «j 1f J 


•a - [A3/(G f - B3>] (7) 

e 12f f 

The inelastic strain rates and the rate of evolution of the resistance 

both in the matrix and fibers are obtained from equations (1) and (2) upon 

using the appropriate material constants. The evaluation of local matrix and 
fiber stress state under strain control requires the simultaneous solution of 

eight differential equations; the first four equations of the set of equa- 

tions (6), equation (7), and two evolution equations of the constituents state 
variables. These equations reduce to five if the fibers are elastic. 


Load Control 

The internal in plane stress state of the composite lamina subjected to 
biaxial stress control shown In figure 2(a) can be obtained by assuming a con- 
tinuum stress field for the composite. In this paper the field is assumed uni- 
form; therefore, the longitudinal, lateral, and shear continuum stress rates 
under the biaxial stress rates shown In figure 2 are given by 

S 11c = S x cos2 0 + S y s1n2 0 
2 ’ 2 

S 22c = S x s1n 9 + S y cos 9 


( 8 ) 


S 12c = - <S x ” S y J s1n 0 cos 0 

where S x and Sy are the applied stress rates in the x and y directions. 
The continuum stresses S-f -| c , S22c and S-j2c are parallel to the continuum 
strains having the same indices. The lateral and shear stress in the fibers 
and matrix are equal to satisfy the equilibrium at fiber surfaces and are 
therefore assumed equal to the continuum lateral and shear stresses at any 
time. The longitudinal fiber stress Is different from that of the matrix. 

The continuum longitudinal stress is an average stress that provides the same 
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longitudinal traction In the composite. Therefore, the longitudinal stresses 
of fibers and matrix are related according to the rule of mixtures, 
s 11m = (1 + v) Sii c - v Sf i f . Notice that the longitudinal fiber and matrix 
total strain rates are equal. Therefore it follows that the fiber longitudi- 
nal stress rate and the stress state in fibers and matrix are given by 


’Ilf 


S 11c * H 2 S 22 c 


H 3 




S 11m = (1 + v) s ii c - v s nf 


S 22m = S 22f ' S 22c 


S 12m " S 12f = S 12c (9) 

where. Hi = Eft(l + v)/(E m + vEf)], H 2 = [(ufE m - u m Ef)/(E m + vEf)], and 
H 3 = [EmEf/(E m + vEf)]. Thus the evaluation of the stress state under stress 
control requires the simultaneous solution of three differential equations. 

The first one of the set of equation (9) and two evolution equations for the 
state variables of the fibers and the matrix. For elastic fibers, the three 
equations reduce to two. 


EXPERIMENTAL PROGRAM 

The material studied was a 4 ply, unidirectional tungsten-fiber reinforced 
copper matrix composite. The volume fraction of the 0.008 in. -diameter tung- 
sten fiber was 9 percent, and the orientation of the fibers was parallel to 
the load axis (i.e., 0 = 0°). Load-controlled fatigue experiments were per- 
formed at 260 °C in vacuum (<1.0E - 5 torr). Five specimens were tested to 
failure and two other tests were Interrupted. 

Figure 3(a) shows the typical cyclic behavior of the composite under 
cyclic tension-tension load control. The composite specimen ratchets continu- 
ously up to failure. Also, the range of cyclic strain decreases continuously 
with cycling, and the hysteresis loops become nearly elastic at failure. Fig- 
ure 3(b) is a plot of the maximum cyclic strain of six specimens as a function 
of the applied number of cycles. The ratcheting behavior Is somewhat analogous 
to typical creep behavior of monolithic materials. The ratcheting rate is high 
at the beginning of the test and decreases to approach a steady state ratchet- 
ing rate after about 50 to 100 cycles. This steady state ratcheting regime 
spans over a quarter of the life for the lowest load level and up to about 
three-quarters of the life for the highest load level. The ratcheting rate 
then increases as failure approaches. The failure strain Increased with 
increasing maximum cyclic stress, ranging from 4.7 percent for a maximum 
stress of 35 ksi to 12.7 percent for a maximum stress of 41.1 ksi. 

Examination of the failed specimens revealed that fatigue cracks nucle- 
ated in the copper matrix via grain boundary cavitation. The cavities linked 
together to form cracks that grew around the fibers. However, some fibers 
broke before the composite specimens failed. No oxidation of the specimens 
was observed. 
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APPLICATION AND DISCUSSION 


The material constants C and D should be evaluated from the fatigue 
data at the temperature of application. Such data were not available for 
tungsten and are not needed for the current application. The copper fatigue 
data at 260 °C were not available either. However, room temperature high 
cycle fatigue data of copper was reported in (ref. 3). The values of C and 
D at 260 °C are assumed to be equal to those at room temperature (table I). 

The initial values of resistance X0 for both copper and tungsten were 
obtained to provide a calculated transient tensile response close to the 
observed one. These values also are given in table I. 

The fatigue shake down resistance of the Insitu composite constituents 
may be different from the resistance of the individual constituents. There- 
fore, the shake down resistance of the matrix is evaluated from the composite 
fatigue data. Specimens 3, 4, and 5 were used to develop a correlation between 
the matrix shake down resistance and applied stress. To develop such correla- 
tion, the tungsten fiber resistance is assumed to be constant and equal to 
its Initial value. The matrix shake down resistance Is obtained Iteratively. 
Initially, the matrix Is assumed linear elastic, and equation (9) is used to 
obtain an initial estimate for L s of the matrix. This estimate is used in 
equation (3) to obtain an Initial estimate for X s of the matrix. Equa- 
tion (9) is then used to update L s , which Is subsequently applied to equa- 
tion (3) to update X s . This process continues until the values of the L s 
and X s for the matrix converge. An effective stress Is defined as the square 
root of the maximum value of 3*32 for the local applied cyclic stress in the 
matrix. Upon correlating the shake down resistance of specimens 3, 4, and 5 
to their effective stress. It can be shown that the copper shake down resist- 
ance at 260 °C is almost constant and equal to 69.0 ksi. This value for shake 
down is applied to predict the lives of the specimens 1, 2, 6, and 7. 

The calculated lives of these specimens are given In table II. Good 
agreement is evident between the observed and predicted low cycle fatigue 
specimens 3, 4, 5, 6, and 7. The agreement Is not good for the high cycle 
fatigue specimens 1 and 2. This can be attributed to the fact that the esti- 
mated value of X s * 69.0 Is obtained from composite low cycle fatigue 
results, which may not be appropriate for high cycle fatigue. 

The variation of predicted life with fiber orientation Is shown In 
figure 4 for the same stresses as had been applied to specimen 7. This predic- 
tion is based on the assumption that composite fatigue failure remains matrix 
controlled and that the matrix fails due to Inelasticity micromechanisms. 
Notice that the effective stress and L s in the matrix decreases as the angle 
of fiber orientation Increases up to an angle of approximately to 33°, after 
which they monotonlcal ly Increase as shown in figure 4. The decrease in L s 
Is followed by a predicted Increase in matrix life. The composite life at any 
given angle of fiber orientation, however, could be fiber-controlled or fiber/ 
matrix interface-controlled. An appropriate life prediction requires the pre- 
dicted lives of each potential contributor on figure 4. The composite life at 
any orientation would be the lowest of the predicted constituent lives. 
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FUTURE RESEARCH 


Three recommendations for future research are proposed. First, the analy- 
sis should be developed to account for the effect of inelasticity due to ther- 
mal strain. Such a step is mandatory If the model is to be applied to predict 
thermomechanical fatigue life. Second, it is required to conduct higher cyclic 
life composite fatigue experiments to Improve the representation of the shake 
down resistance for high cycle fatigue. Third, composite fatigue experiments 
should be conducted at different angles of orientation to verify the validity 
of the proposed life prediction method. 
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TABLE I. - MATERIAL CONSTRUCTION 


Cons 

Uni ts 

Cu 

W 

AO 

1/sec ksi 

10 000 

10 000 

r 


7.0 

15.00 

C 

1/ksi 

-.2526 


D 


-37.2 


xO 

ksi 

23.0 

800.0 


TABLE II. - FATIGUE RESULTS 
[6 - 0 °] 


Number 

Maximum 

stress, 

ksi 

Minimum 
stress , 
ksi 

Nf. 

exp 

N f , 

theo 

1 

33.5 

1.42 

a 30 419 

28 430 

2 

34.4 

1.53 

*62 403 

18 510 

3 

35.0 

1.48 

b 12 588 

12 370 

4 

37.2 

1.53 

b 3 117 

4 648 

5 

38.6 

1.52 

b 2 727 

2 464 

6 

40.2 

1.49 

1 864 

1 250 

7 

41.1 

1.51 

933 

8 855 


a Test interrupted prior to failure 
“Test used to calculate material constants 
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LOADING FUNCTION, L RESISTANCE TO FLOW, 


x 



FIGURE 1. - SCHEMATIC VARATIONS IN THE RESISTANCE OF 
MATERIAL TO FLOW AND IN LOADING FUNCTION WITH THE 
NUMBER OF APPLIED FATIGUE CYCLES. 



(a) APPLIED BIAXIAL STRESS AND (b) FIBER AND MATRIX IN PLANS 
STRAIN RATES. STRESSES AND STRAINS. 


X CONTINUUM DISPLACEMENT 

v 



(C) FIBER, MATRIX, AND CONTINUUM DISPLACEMENTS AND STRAINS IN 
A DIRECTION PERPENDICULAR TO FIBERS. 

FIGURE 2. - ELEMENT OF COMPOSITE LAMINA. 


MAX STRESS = 41.10 Ksi 
MIN STRESS = 1.51 Ksi 


(a) VARIATION OF MATRIX EFFECTIVE STRESS WITH FIBER 
ANGLE. 
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